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ABSTRACT 

We use all available fields with deep NICMOS imaging to search for Jhq- 
dropouts (i/i6o,As ^ 28) at z ~ 10. Our primary data set for this search were 
the two Jiio + Ri 60 NICMOS parallel fields taken with the Advanced Camera for 
Surveys (ACS) Hubble Ultra Deep Field (HUDF). The 5a limiting magnitudes 
were 28.6 in Jhq and 28.5 in Hiqq (0.6" apertures). Several shallower fields were 
also used: Jno + Rieo NICMOS frames available over the Hubble Deep Field 
(HDF) North, the HDF South NICMOS parallel, and the ACS HUDF (with 5 ct 
limiting magnitudes in Jno and Rieo ranging from 27.0 to 28.2). The primary 
selection criterion was {Jiio—Hiqq)ab > 1-8. Eleven such sources were found in all 
search fields using this criterion. Eight of these were clearly ruled out as credible 
z ^ 10 sources, either as a result of detections (> 2a) blueward of Jno or their 
colors redward of the break {Hiqq — K ~ 1.5) (redder than >98% of lower redshift 
dropouts). The nature of the 3 remaining sources could not be determined from 
the data. The number appears consistent with the expected contamination from 
low-redshift interlopers. Analysis of the stacked images for the 3 candidates 
also suggests some contamination. Regardless of their true redshifts, the actual 
number of 2 ; ~ 10 sources must be < 3. To assess the significance of these results, 
two lower redshift samples (a z ~ 3.8 R-dropout and 2 : ~ 6 i-dropout sample) 
were projected to ^ ~ 8 — 12 using a (1 + z)~^ size scaling (for fixed luminosity). 

They were added to the image frames, and the selection repeated, giving 15.6 
and 4.8 Jno-dropouts, respectively. This suggests that to the limit of this probe 
(«i 0.3L*^3) there has been evolution from 2 ; ~ 3.8 and possibly from z ^ 6. This 
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is consistent with the strong evolution already noted at z ^ 6 and ~ 7.5 relative 
to z ^ 3 — 4. Even assuming that 3 sources from this probe are at z zz 10, the 
rest-frame continuum UV (~ 1500A) luminosity density at ~ 10 (integrated 
down to O.SL^^g) is just 0.19lg;Q9X that at z ^ 3.8 (or 0.19 ](q;JqX including the 
small effect from cosmic variance). However, if none of our sources is at ~ 10, 
this ratio has a la upper limit of 0.07. 

Subject headings: galaxies: evolution — galaxies: high-redshift 


1. Introduction 

The detection of galaxies at the highest redshifts {z > 7) continues to be a difficult 
endeavor. Due to the redshifting of UV light into the infrared (IR) and the well-known 
limitations of current IR instruments, searches for these objects require almost prohibitive 
amounts of telescope time. Nevertheless, small amounts of deep IR data do exist, and they 
can be used to set constraints on very high redshift galaxies. One notable example is the 
2 ; 85 o-dropout sample compiled by Bouwens et al. (2004c) which, although limited by small 
numbers (~ 5 objects) and £eld-to-£eld variations, provided a £rst detection of galaxies at 
redshifts beyond z ^ 7. In this paper, we look at the prevalence of galaxies at 8-12 by 
applying the dropout technique to the wide variety of deep FllOW and F160W-band helds 
that have been imaged with NICMOS (hereinafter, refered to as Jno and Hiqq, respectively). 
Our principal data set in this study will be the two deep NICMOS parallels taken with the 
UDF (each has ~ 200 orbits of data), but we will complement this held with a variety of 
shallower fields possessing similar Jug -|- Tfieo imaging. These fields include the HDF-North 
Thompson field (Thompson et al. 1999), the HDF-North Dickinson field (Dickinson et al. 
1999), the deep HDF-South parallel NICMOS field (Williams et al. 2000), and the NICMOS 
footprint on the UDF itself (Thompson et al. 2005). Galaxies in the redshift range z ^ 3 —12 
are currently of great interest due to indications that reionization of the universe may have 
started as soon as ~ 17 ± 5 (Kogut et al. 2003) and that galaxies may have played a 
major role in this process (e.g., Yan & Windhorst 2004; Stiavelli et al. 2004). We take 
to denote the characteristic luminosity of galaxies at z = 3 (Steidel et al. 1999). AB 
magnitudes are used throughout. We assume (Dm, Da, h) = (0.3, 0.7, 0.7) (Bennett al. 2003). 
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2. Observations 

(a) NICMOS parallels to the UDF. The two NICMOS parallels to the UDF (taken 
during the first and second epochs of ACS observations) make up our primary data set. 
Each parallel consisted of ~200 orbits of data, split between two pointings which overlap 
by ~25% of a NIC3 pointing (~ 3.0 x 10^ s were taken at the hrst and ~ 8.6 x 10"^ s at 
the second). The observing time at each position was split nearly equally between Jno 
and hfi6o-band observations. While reductions of these helds are available through STScI 
(http://archive.stsci.edu/prepds/udf/udf_hlsp.html), we used the NICMOS procedures de¬ 
scribed in Thompson et al. (2005) to perform our own reduction. The deeper of the two 
pointings making up each parallel had approximate 5a depths of 28.6 and 28.5 (0.6" aperture) 
in the Jno and i^igo-bands, respectively. The shallower pointings were some ~0.6 mags less 
deep. Our reductions were drizzled onto a 0.09"-pixel scale, with approximate PSFs of 0.32" 
and 0.34" FWHM in the Jno and Ffigo-bands, respectively. For optical coverage on these 
helds, the Woe and 2 : 85 o-band images from Galaxy Evolution Morphology Survey (Rix et al. 
2004: images #25, #42, and #49) were used (Blakeslee et al. 2003). 

(b) Shallower NICMOS Fields. To add area at brighter magnitudes, we included a 
number of other helds with deep Jno -l- Jfieo imaging in our Jno-dropout search. A list of 
these helds is provided in Table 1, together with their approximate 5a limiting magnitudes 
(0.6" aperture) and selection areas. Other multi-wavelength data for these helds include 
deep UBVI coverage for the HDF-N (Williams et al. 1996), ultra deep ACS BViz coverage 
for the UDF (Beckwith et al. 2005, in preparation), NICMOS A '222 coverage of the HDF- 
S NICMOS parallel (Williams et al. 2000), and deep IRAC data in the 3.6/im and 4.5yum 
channels for the GOODS North and South (Dickinson et al. 2005, in preparation). 


3. Analysis 

SExtractor (Bertin & Arnouts 1996) was used to do object detection and photometry. 
Detection was performed using an aggressive 2a detection threshold on the iJieo-band images. 
Our catalogs were then cleaned of contamination from more extended background artifacts by 
requiring objects to be 5a detections in a 0.5"-aperture. We used scaled aperture magnitudes 
(Kron 1980: SExtractor MAG_AUTO) for our Jfigo-band total magnitudes. Colors were also 
measured with scaled apertures, but with a much smaller Kron factor to maximize the S/N. 
For images with diherent PSFs, colors were only measured after the higher resolution image 
had been degraded to match the broader PSF and only in a circular aperture which was at 
least 2x the FWHM of the broader PSF (see Bouwens et al. 2003 and B05 for more details). 



(a) JiiQ-dropout selection. Our principal selection criterion for our sample was a (Jno — 
Hiqo)ab > 1-8) color cut, where the Jno-band flux was replaced by its la upper limit in cases 
of a Jiio-band non-detection. This criterion was chosen to minimize contamination from all 
objects except z ^ 2 — 5 evolved galaxies (Figure 1). In addition, objects were required to 
be non-detections (< 2a) in all optical bands as a result of absorption from the intergalactic 
medium and to be blue (< 1 mag) in all Hiqq — K ot Hiqq — I RAC colors, since only a small 
fraction (< 2%) of star-forming galaxies have sufficiently red UV continuum slopes f3 (i.e., 
f3 > 0.5) to produce these colors (Adelberger & Steidel 2000; Schiminovich et ah 2004). 

Applying our (Jno — Hiqo)ab > 1-8 color cut to object catalogs from both our primary 
and shallow search helds, 11 objects were found (Table 2). Of the 11, 6 were readily detected 
(> 2a) in the bluer optical bands. Two had very red Hiqq — K colors (~ 1.5) and so 
also appear to be low-redshift interlopers. This was a little uncertain for one of the two 
(HDFSPAR-48278437) due to the marginal nature of its iF 222 -band detection, but IRAC 
imaging should clarify this issue. The hnal three sources (shown in Figure 2 with two red 
objects from our selection) could not be excluded by either criterion, though this may be 
largely due to the fact they were in the UDF-parallel helds. The optical imaging in these 
helds are not particularly deep, nor is there any deep Spitzer or K-band data to measure 
colors longward of the break. 

(h) Low-Redshift Contamination. To help determine whether the 3 potential 2 : ~ 10 can¬ 
didates were low-redshift interlopers or not, we made simple estimates of the contamination 
from diherent types of sources. The only known stellar sources red enough to match these 
objects are extreme carbon stars or Mira variables (Whitelock et ah 1995) but contamina¬ 
tion from such sources seem unlikely due to both their rarety and high intrinsic luminosities, 
which would put them well outside the Galaxy (Dickinson et al. 2000). Contamination from 
redder, evolved galaxies is difficult to estimate given the large uncertainties on the shape 
of the LF of these galaxies at ~ 2 — 3. Here, we ignore these complications and simply 
scatter the colors of a brighter set of galaxies (iPieo.As ~ 24 — 26) from our helds to fainter 
magnitudes. Performing these experiments on the faint source population from the deep 
parallels, we estimate ~ 1 — 3 such contaminants. This seems consistent albeit a little less 
than the ~ 3 sources obtained. A stack of the Vqqq and Jno exposures for all 3 sources 
showed detections of 0.8(j and 1.2a, respectively, again snggesting some contamination. 

(c) Expected Numbers/Incompleteness. Having fonnd at most 3 possible 2 : ~ 10 sonrces 
over our entire 14.7 arcmin“^ search area, it is interesting to ask how many we might have 
expected assuming no-evolution from lower redshift. Two different redshift samples are 
considered as baselines: (1) a 2 : ~ 3.8 R-dropout sample from GOODS (B05) and (2) a z ^ 6 
i-dropout sample from the HUDF (Bonwens et al. 2005b). As in other recent work (Bonwens 
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et al. 2004a,b,c), we project this sample to higher redshift (over the range z ^ 7 — 13: see 
Figure 3) using our well-established cloning machinery (Bouwens et al. 1998a,b; Bouwens 
et ah 2003; BOS), accounting for pixel morphologies, the angular-size distance relationship, 
NICMOS PSFs, k-corrections, and object-by-object volume densities. Transformed objects 
are added to the present NICMOS data and the selection procedure repeated. Direct use of 
the data appears to be the best way of accounting for the substantial variations in S/N which 
occur across NICMOS mosaics while including possible blending with foreground galaxies. 
Simulations were run over ~30x the area of the helds. Assuming no-evolution in size, 5.7 
and 4.1 objects are expected in total (over the 6 helds) for our 2 : ~ 3.8 B and z ^ 6 Fdropout 
samples, respectively. If we account for the increase in surface brightness expected due to 
the ~ (1 -|- z)~^ size scaling observed at 2 < 2 : < 6 (Bouwens et al. 2004a,b; Ferguson et 
al. 2004), the expectations increase to 15.6 and 4.8, respectively. Steeper size scalings (e.g., 
(1 - 1 - z)~'^) yield 21.2 and 6.3 Jno-dropouts, respectively, while use of a bluer UV continuum 
slope f3 (e.g., f3 ~ —2.4) resulted in 16.2 and 4.0, respectively (assuming a (1 -|- z)~^ size 
scaling). We note that the two NICMOS parallels to the UDF account for ~ 73% of the 
expected numbers (and are therefore our primary constraints), though each of the 6 helds 
contributes at least a few percent. We take the ( 1 - 1 - 2 ;)“^ scaling estimates, 15.6 from 2 ; ~ 3.8 
and 4.8 from z ^ 6 as the most likely expected values. 

(d) Previous Work. One previous study of > 8 candidates was carried out by Yahata 
et al. (2001) on the 0.8 arcmin^ NIC3 parallel to the HDF South featured here. In that 
paper, 8 2 : > 10 candidates were reported (5 of which were selected in the Tfigo-band) based 
upon 9-band (UBVRIJHK-I-STIS) photometric redshifts. What became of these 5 candidates 
here? Looking through our catalogs for this held, we found that one (SB-NI-0915-0620) had 
Jiio - Hiqo colors ( 0.4) that were clearly inconsistent with a high redshift identihcation. 

The other 4 appeared to be too faint {Hiqq ab ^ 27.5) to set strong lower limits on the 
(■^iio — Him)AB colors and thus make robust statements about their redshifts. 


4. Implications 

We have carried out a search for z ^ 10 J-dropouts and found at most 3 possible 
candidates. Since 15.6 ±3.9 and 4.8 ±2.2 candidates are expected (Poisson errors) based on 
a (1 ± z)~^ scaling of 2 : ~ 3.8 and 2 ~ 6 samples, this suggests there has been appreciable 
evolution at the bright end of the luminosity function. Since an ±*=3 galaxy at 2 ~ 10 has 
a i7i6o,AB-band magnitude of 26.7, our search (hfieo.As ^ 28) tells us something about the 
galaxy luminosities brightward of 0 . 3 L*^ 3 . For simple Poissonian statistics and assuming 
that all 3 sources are at 2 ~ 10, the current hndings are inconsistent with no-evolution at 
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the 99.98% and 71% confidence levels, respectively (equivalent to 3.8(j and la). 

Of course, cosmic variance is bound to be important for fields of this size. Assuming 
a selection window of width Az = 2.5, a ACDM power spectrum, and a bias of 7-which 
corresponds to the rough volume density of sources ~ 10“"^ Mpc“^ explored by this probe 
(Mo & White 1996)-we calculate ~27% RMS variations field-to-field for single 0.8 arcmin^ 
NIC3 pointings (using a pencil beam window function). Since our deepest two fields provide 
the primary constraints and they are essentially independent, the total number of J-dropouts 
found here is expected vary by ~19% RMS relative to the cosmic average. Thus, we expect 
~16 ± 5 and ~5 ± 2 dropouts, respectively, for simple projections of our lower redshift 
samples, and so the present findings are inconsistent with no-evolution at the 99.9% and 
68% confidence levels, respectively (equivalent to 3.3a and 1.0a). 

It is conventional to cast these findings in terms of the rest-frame continuum UV (~ 
1500A) luminosity density. We first estimate the luminosity density under the assumption 
that the 3 candidates (Figure 2) are at z ^ 10. The result is puv(z 10)/puviz ~ 3.8) = 
O.lOt^og (3 over 15.6 for Poissonian statistics) or puv{z ~ ^0)/puv{z ~ 3.8) = O.lOj^o^Q 
(including cosmic variance). Taking the value of the luminosity density at 2 : ~ 3.8 (Giavalisco 
et al. 2004) integrated down to 0.3L*^3, we find a UV luminosity density puv{z 10) = 
2.5l};2 ^ Hz“^ (Poissonian statistics) or puv{.z ~ 10) = x 10^^ erg s“^ Hz“^ 

(including cosmic variance). However, it is quite plausible that none of these candidates are 
at z zz 10. The expected high contamination level combined with the now well established 
changes in the LF between 2 ; ~ 6 — 7 and z 3 suggests that it is more likely that these 
sources are low redshift objects. Assuming that none are at 2 : ~ 10, the la upper limit 
is: puv{z ~ 10)/puv{z ~ 3.8) < 0.07 (both for simple Poissonian statistics and including 
cosmic variance). In terms of the luminosity density, this limit is puv{z ~ 10) < 0.9 x 10^® 
erg s“^ Hz“^. Adopting a Salpeter IMF and using canonical relations to convert this into 
a star formation rate density (Madau et al. 1998), we can plot the present determination 
(integrated down to 0.3L*^3) against previous determinations at lower redshift (Figure 3). 
The observations now more clearly than before allow us to trace the number of f/U-bright 
systems over the interval 0 < 2 < 10. The space density of high luminosity systems seems 
to peak at 2 ~ 2 — 4 and decline fairly rapidly to both higher and lower redshift. 

As we conclude, it is somewhat sobering to realize that ~740 orbits of deep NICMOS 
imaging went into this search and only 3 possible 2 ~ 10 candidates were found (all of which 
may be at low redshift), showing how difficult it is to map out the formation of galaxies 
at these early times with current technology. It is exciting nevertheless to realize that in 
the future these surveys will be executed much more efficiently. For example, surpassing the 
current NICMOS data set would require just ~23 orbits with HST WFC3 and ~ 1000 — 2000 



seconds with JWST. 
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Table 1. 


Jiio-dropout Search Fields. 



Area 

5(7 limit 

Field 

(□') 

Rio 

Hieo 

HDF-N Thompson 

0.8 

27.8 

28.1 

HDF-N Dickinson 

5.2 

27.2 

27.0 

HDF-S NIC-Par 

0.8 

28.2 

28.2 

UDF NIC-Par-1 

1.3 

28.6 

28.5 

UDF NIC-Par-2 

1.3 

28.6 

28.5 

UDF Thompson 

5.5 

27.7 

27.5 


Table 2. 


Red (Jiio - Hiqq)ab > 1-8 objects 


in our search helds.®" 


Object ID 

R.A. 

Decl. 

Hiqo 

J - H 

V - H 

z- H 

H - K 


rw(") 

UDF-39188323 

03:32:39.18 

-27:48:32.3 

23.3±0.1 

2.0 

3.8 

2.5 

0.4 

1.2 

0.26 

UDF-42888094 

03:32:42.88 

-27:48:09.5 

24.8±0.1 

2.7 

4.3 

3.3 

1.0 

2.9 

0.25 

UDFNICPARl-01191115 

03:33:01.23 

-27:41:11.5 

27.0±0.1 

1.8 

0.9 

>-0.2 

— 

— 

0.24 

UDFNICPARl-04151142 

03:33:04.18 

-27:41:14.1 

28.2±0.2 

1.8 

>-0.0 

>-1.3 

— 

— 

0.18 

UDFNICPARl-05761077 

03:33:05.80 

-27:41:07.6 

27.8±0.2 

>2.0 

>0.3 

>-1.0 

— 

— 

0.21 

UDFNICPAR2-07352493 

03:33:07.35 

-27:52:49.3 

27.9±0.2 

1.8 

0.6 

>-0.7 

— 

— 

0.14 

UDFNICPAR2-09593048 

03:33:09.59 

-27:53:04.8 

27.3±0.2 

1.9 

>0.9 

>-0.4 

— 

— 

0.19 

HDFN-46242594 

12:36:46.24 

62:12:59.4 

26.5±0.1 

>2.1 

1.1 

— 

<1.7 

— 

0.17 

HDFN-51763567 

12:36:51.76 

62:13:56.7 

25.0±0.1 

2.0 

2.1 

— 

0.6 

— 

0.31 

HDFN-56253227 

12:36:56.25 

62:13:22.7 

25.2±0.1 

>3.2 

>4.2 

— 

1.5 

2.0 

0.20 

HDFSPAR-48278437 

22:32:48.27 

-60:38:43.7 

26.0±0.1 

>2.5 

— 

— 

1.5 

— 

0.32 


^All astrometry uses the J2000 equinox. Magnitudes are in the AB system. All limits are 2a except the limits on the Jiio — 
color which are only la. Total magnitudes and colors are derived in Kron apertures of different sizes (B05). Object IDs are the last four 
significant figures given in this table for the Right Ascension and Declination. 
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Fig. 1.— {JiiQ — Hiqq)ab color versus redshift for a number of different SEDs. Shown are the 
Coleman, Wu, & Weedman (1980) elliptical template {red line), Sbc template {blue line), 
Scd template {green line), and different reddenings {E{B — V) = 0.0,0.15,0.3) applied to 
a 10® yr burst {black solid, dotted, and dashed lines, respectively). Extremely red objects 
such as those found in recent IRAC selections (e.g., Yan et ah 2004) have very similar colors 
(in the infrared) to 2.5 Gyr bursts and are included here as the dashed red line. The red 
arrow near z ~ 0 denotes the range of colors expected for low mass stars (Knapp et al. 
2004). Our (Jiio — Hiqq)ab > 1-8 Jno-dropout criterion is shown as a dotted horizontal 
line and provides a good balance between minimal contamination (selecting evolved SEDs 
at intermediate redshifts z ~ 2 — 5) and selecting objects at z ^ 10. 
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V J H 

30.2+0.8 >28.8 25.6+0.1 

il 

HDFN-56253227 
>29.3 26.9+0.1 

HDFSPAR-48278437 


>29.1 

>29.5 27.7+0.2 


A.. A 


UDFNICPAR2-09593048 

W >29.1 

■^30.1 ' ^^1+0.2 


r 

UDFNICPARl-057^077 



29.1+0.9 


UDFNicPARl-0415Tl42 


Other Bands 
23.6+0.1 



IRAC 3.6!im 


25.2+0.5 


NICMOS K 


Fig. 2.— Postage stamps images (Jno^^ieo bands) of 5 objects which met our (Jno — 
Hiqq)ab > 1-8 selection. The leftmost and rightmost columns show postage images in 
a bluer (Vgoe band) and redder band (the IRAC 3.6/rm channel for HDFN-56253227 and 
NICMOS K222 for HDFSPAR-48278437), respectively. The overplotted magnitudes (and la 
upper limits) were measured within a 0.6''-diameter aperture. The object in the top row 
(HDFN-56253227) is the well-known Dickinson et al. (2000) Jno-dropout. Its extremely red 
colors ((iPieo — K)ab = 1-5 and (iPieo — I ab = 2) suggest that it is at low redshift 
( 2 : ~ 2 — 3). The second object (HDFSPAR-48278437) was similarly excluded from our list 
of 2 : ~ 10 candidates because of its marginal detection (si 2 a) in the NICMOS iF 222 -band 
and (i7i6o — K 22 o)ab color of 1.5. The hnal three objects could be at 2 : ~ 10, but await 
limits at longer wavelengths (from Spitzer) and good optical data to exclude the possibility 
that they are highly reddened or evolved galaxies. Objects UDFNICPAR2-09593048 and 
UDFNICPARl-04151142 appear to be just detected in the Jno-band at the la level, so this 
may indicate they are low redshift interlopers. The hnal object (UDFNICPARl-04151142) 
is just marginally resolved, but the only stellar sources red enough to match its colors would 
appear to be too bright (§3b). Each postage stamp is 2.9"x2.9" in size. 







13 


t (Gyr) 



26 

25.5 
25 

24.5 


o 

OQ 

ts 

5 

o 

o 


(B 

-i 

OQ 

a> 

m 


I 

N 


■ra 

o 



Fig. 3.— The cosmic star formation rate density versns redshift with no extinction correction. 
The star formation rate density (integrated down to 0.3 - the limit of onr ^ ~ 10 search) 

was calcnlated from the Inminosity density in the rest-frame UV continnnm (~ ISOOA: 
plotted on the right vertical axis) using canonical assumptions (e.g., Madau et ah 1998) 
and a Salpeter IMF. The open red square at z ~ 10 shows our result if 3 objects from 
this study prove to be at z ~ 10, while the large downward pointing arrow shows our Icr 
limits if none are (see text). Included on this plot are also estimates by Schiminovich et 
ah (2004) (open black squares), Steidel et ah (1999) (green crosses), Giavalisco et ah (2004) 
(black diamonds). Bunker et ah (2004) (solid blue square), Bouwens et ah (2004c) (solid red 
square), and Bouwens et al. (2005b) (solid red circle). Consistent with past practice, the 
error bars reflect poissonian uncertainties. Large scale structure (cosmic variance) would 
add an estimated ±20% for many of the lower redshift points (e.g., Somerville et al. 2004), 
±50% for the ~ 7.5 point, and ±19% for the ~ 10 point. Lower left: The surface density 
of Jiio-dropouts predicted from a (1 ± z)~^ size scaling of a 2 ; ~ 3.8 i?-dropout sample from 
GOODS (B05) at the depths of the NICMOS parallels to the UDF (shaded red regions). The 
blue shaded regions show the equivalent surface density for all our search helds. The drop 
in the predicted counts at fainter magnitudes arises from incompleteness. Lower right: The 
redshift distribution for objects satisfying our Jno-dropout criterion in the above simulations 
(after distributing our “cloned” objects over the interval 2 ; ~ 7 — 13 in proportion to the 
cosmological volume element). The primary conclusion to be drawn from these results is 










